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ABSTRACT 


An  hi  radio  experiment  was  performed  to  measure  the  high-altitude 
effect  of  the  vertically  traveling  pressure  wave  resulting  from  a  large 
ground-level  explosion.  The  blast--Pro ject  Snowball--consisted  of  500  tons 
of  TNT  and  was  detonated  at  the  Suffield  Experimental  Station,  Alberta, 
Canada,  July,  1964.  The  ionospheric  disturbance  was  monitored  using 
vertical-incidence,  phase-sensitive  sounders  located  85  km  from  ground 


Simple,  linear,  acoustic  theory  was  used  to  calculate  the  onset  time, 
amplitude,  and  period  of  the  radio-signal  disturbance.  These  calculations 
agree  closely  with  measurements  taken  by  the  vertical-incidence  sounders-- 
onset  time  was  predicted  within  10  sec,  and  both  amplitude  and  period 
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agreed  within  a  factor  of  two. 
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I .  INTRODUCTION 


On  17  July  1964  at  1058  MST ,  approximately  500  tons  of  TNT  was 
detonated  on  the  ground  surface  at  the  Suf field  Experimental  Station 
(SES),  Ralston,  Alberta,  Canada  (Fig,  l)’.  The  detonation  culminated  a 
comprehensive  program  of  blast  studies,  under  the  sponsorship  of  the 
Canadian  Defence  Research  Board,  under  a  tripartite  arrangement  with 
Canada,  the  United  Kingdom,  and  the  United  States  represented.  The  U.  S. 
participation  in  this  test  was  named  Operation  Snowball.  Stanford's 
primary  objective  in  participating  in  this  program  was  to  study  the 
effect  of  the  blast  on  the  ionosphere. 

The  shock  wave  initially  generated  by  the  atmospheric  explosion 
degenerates  rapidly  into  a  low-frequency  sound  wave,  which  then 
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FIG.  1.  LOCATION  OF  SNOWBALL  EXPLOSION. 
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propagates  to  great  distances  with  little  attenuation.  Acoustic  energy 
directed  vertically,  for  example,  can  reach  the  ionized  regions  of  the 
extreme  upper  atmosphere.  The  effect  of  such  pressure  disturbance  on 
the  ionosphe  ’e  may  be  studied  from  the  ground  using  high-frequency  (hf ) 
radio  signals  that  reflect  vertically  from  the  ionized  layer  directly 
above  the  blast  (Fig.  2).  The  disturbance  causes  localized  electron- 
density  variations,  which,  in  turn,  influence  the  phase  path  of  the 
reflecting  radio  wave. 

For  the  purpose  of  instrumenting  the  ionosphere  above  the  Snowball 
explosion,  two  vertical -incidence ,  phase-sensitive  ionospheric  sounders 
were  located  50  miles  magnetic  north  and  south  of  the  blast — ne^r  the 
towns  of  Arneson  and  Bow  Island  (Fig.  l) —  in  an  attempt  to  measure  the 
influence  of  the  earth's  magnetic  field  on  acoustic-wave  propagation. 
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FIG.  2.  CROSS  SECTION  SHOWING  EXPLOSION  AND  GEOMETRY  OF 
VERTICAL-INCIDENCE  PULSE  SOUNDER. 
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I I .  THEORY 


The  following  section  presents  a  brief  outline  of  the  method  ured 
to  predict  the  disturbance  detected  by  ionospheric  sounders  operating 
no:?*  the  500-ton  Snowball  explosion.  The  theory  is  presented  in  two 
oarts:  a  calculation  of  the  size  and  shape  of  the  pressure  wave  pro¬ 

duced  in  the  ionosphere  by  the  blast,  and  a  prediction  of  the  effect  of 
such  a  pressure  disturbance  on  an  hf  radio-wave  reflection  from  the 
ionosphere . 

A.  PROPAGATION  OF  EXPLOSION-INDUCED  SHOCK  WAVES 

The  mathematical  description  of  the  propagation  of  acoustic  waves 
generated  by  explosions  in  the  atmosphere  is  facilitated  by  considering 
separately  near-  and  far-field  effects.  The  shock  wave  traveling  in 
the  near  field  must  be  described  by  a  set  of  nonlinear  differential 
equations.  For  the  purposes  of  this  paper,  the  blast-wave  computations 
of  Brode  [Ref.  l]  were  used  to  describe  near-field  effects.  Brode ' s 
results  extend  from  the  surface  of  the  explosive  out  to  a  point  at  which 
the  maximum  shock  overpressure  has  been  reduced  to  6  percent.  Figure  3 
shows  the  6-percent  wave  shape  given  for  a  500-ton  blast. 


34268 

FIG.  3.  PRESSURE  VS  DISTANCE  FROM  SOURCE  FOR  SHOCK  WAVE  AT 
6 -PERCENT  OVERPRESSURE. 
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From  this  point  outward,  the  disturbance  may  be  treated  as  a  linear 
acoustic  wave.  To  simplify  calculations,  an  approximation  to  Brode's 
waveform  (Fig.  4)  was  assumed  uniformly  distributed  over  a  hemispherical 
wavefront  1150  meters  from  the  source  (Fig.  5).  The  following  discussion 
summarizes  the  linear  acoustic  theory  used  to  describe  propagation  in 
the  far-field  regions. 


0.08 
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FIG.  4.  STRAIGHT-LINE  APPROXIT  AT ION  OF  6-PERCENT  OVERPRESSURE  USED 
IN  CALCULATIONS  <  >F  THIS  REPORT. 


1150m - H 


34270 

FIG.  5.  6-PERCENT  WAVEFRONT  USED  IN  CALCULATIONS  OF  THIS  REPORT 
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The  maximum  power  density  carried  by  a  linear,  plane,  harmonic 
acoustic  wave  in  the  atmosphere  is  given  the  equation  [Ref.  2] 


P 


/P 

[  max 

V  p, 


(1) 


where 


p  =  atmospheric  density 

a  =  velocity  of  sound 

7  =  ratio  of  specific  heats,  cp/cv 

p  =  maximum  pressure  in  the  wave 
max 

p  =  ambient  atmospheric  pressure, 
o 


One  may  relate  the  acoustic  properties  of  a  wave  propagating  from 
point  1  to  point  2  in  the  atmosphere  in  terms  of  the  properties  at 
each  point  by  the  simple  ratio 


©'  ■  0  mi  %  • 


(2) 


where  it  is  the  peak  normalized  overpressure,  (p  -  p  )/p  . 

max  o'  o 

This  equation  is  the  basis  of  the  linear  acoustic  theory  used  in 

this  paper.  The  wave  at  point  1  is  taken  as  the  6-percent  wavefront 

and  associated  wave  shape  described  by  Erode.  The  wave  properties  at 
point  2--the  ionosphere  directly  above  the  blast--are  then  calculated 
by  using  Eq.  (2),  Inspection  of  the  equation  shows  immediately  that  all 
terms  are  defined  except  the  maximum  power-density  ratio,  Pg/P^ ' 
values  of  p,  a,  and  7  have  been  taken  from  U.  S.  Standard 
Atmosphere  [Ref.  3].) 

As  a  shock  wave  propagates  upward  in  the  atmosphere,  the  power 

density  in  the  wave  will  be  decreased  by  two  separate  mechanisms,  each 

of  which  may  be  considered  independently.  The  first  loss  is  caused  by 
spreading.  Since  the  wavefront  area  becomes  increasingly  large  with 
propagation  away  from  the  source,  the  total  energy  in  the  wave  is  spread 
over  a  larger  area  and  both  energy  density  and  power  density  are  reduced. 
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The  second  loss  is  introduced  by  atmospheric  attenuation;  in  fact,  the 
atmosphere  acts  as  a  lossy  bandpass  filter  to  sound  waves. 

Taking  these  effects  into  account,  the  power  loss  may  be  expressed 
as 


P 

P 


2 

1 


=r.  AP 

s 


(3) 


where 


APg  =  power-density  loss  factor  caused  by  spreading 

AP^  =  power-density  loss  factor  caused  by  atmospheric 
filtering . 


Shock  waves  generated  by  explosions  in  a  uniform  atmosphere  are 

spherical  and  remain  spherical  as  the  shock  propagates  outward;  the 

power  density  is  inversely  proportional  to  the  area  of  the  front.  Thus, 

the  spreading  loss  factor  APg  (Eq.  3),  for  a  disturbance  propagating 

from  the  3rode  6-percent  radius  (i.e.,  1150  meters  from  the  source)  to 

the  upper  ionosphere  (i.e.,  F-layer  at  220-km  altitude)  in  a 

2  —6 

homogeneous  atmosphere  is  AP^  =  (1.15/220)  =  5.23  X  10 

For  the  real  atmosphere,  the  increased  loss  caused  by  nonuniform 
spreading  was  estimated  with  the  use  of  acoustic  raytracing  (Fig.  6) 
and  found  to  be  about  0.3  times  the  value  given  for  a  homogeneous 
atmosphere.  Thus 


AP 

s 


1.57  X  10 


The  atmospheric-filtering  loss  factor  AP^  (Eq.  3)  is  obtained  from 
the  general  wave  equation  for  linear,  plane  acoustic  waves. 


- -  (x,t)  =  -ma-X — -exp[-ax]exp  j2nf  |t-(x/a )} 

p„  P„  L  *  *J 


(5) 


where  a  =  attenuation  constant  (nepers/meter) 
f  =  acoustic  frequency 
a  =  speed  of  sound 
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FIG.  6.  ACOUSTIC  RAYTRACING  FOR  GROUND-LEVEL  EXPLOSION— RAYS  AT 
10 -DEG  INCREMENTS. 


The  total  power  attenuation  (in  db)  as  a  function  of  frequency  for  a 
wave  traveling  between  two  points  x  and  x  is  then 


AP(f)  =  8.68 


dx 


db 


(6) 


The  attenuation  constant  a  is  a  function  of  frequency  and  altitude 
and  was  shown  by  Rayleigh  [Ref.  2]  to  be 


a(f,h) 


where  T)  =  kinematic  viscosity 

k  =  thermal  conductivity 

7  =  ratio  of  specific  heats,  C  /c 

p  v  ' 


(7) 


7 


SEL-65-062 


Ignoring  the  generally  much-smaller  conductivity  term  and  substituting 
[Ref.  2] 


i)  =  0.499 


v2  1.497  a2 


where  v  -  mean  particle  velocity 
v  =  viscosity  coefficient, 


yields 


a(f,h)  - 


2  2 

2  4n  f 

}  va 


Substituting  this  value  in  Eq .  (6)  gives 


Ap(f)  =  3.42  X  10 2  f  — 
J  va 


where 


AP(f)  =  Kx_2  f" 


K  =  3.42  X  10 


J  va 


The  value  of  K  was  evaluated  graphically  with  the  point  x 
1  1 

taken  as  the  top  of  the  6-percent  wavefront  1150  meters  above  the 
source  and  x2  as  the  F  layer  at  220-km  altitude.  The  resulting  value 


K  -  6.01  X  103  sec"2 

l  “ 


Equation  (9)  above  describes  the  attenuation  of  high  acoustic  fre¬ 
quencies.  Propagation  of  low  frequencies  is  limited  by  the  cutoff 
frequency  f  ,  which  is  caused  by  the  variation  in  atmospheric  density 
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with  altitude.  The  value  of  the  lower  limit  may  be  obtained  by  observing 
that,  to  a  vertically  traveling  sound  wave,  the  atmosphere  has  the 
appearance  of  an  exponentially  tapered  transmission  line.  The  acoustic 
impedance  is  [Ref.  2] 


z  =  ap  ,  (12) 

where  a  is  the  velocity  of  sound  and  p  is  the  atmospheric  density. 
Although  both  the  speed  of  sound  and  density  vary  with  altitude,  varia¬ 
tions  in  the  latter  are  much  more  pronounced,  being  approximately 
exponential.  Thus, 

p  =  p  exp[-(h-h  )/h] 
o  o 

=  p^  expfh^/H]  exp[-h/H]  (13) 

=  p  exp[-h/H]  , 

where  h  =  altitude 

H  =  scale  height 

p  =  density  at  some  altitude  h  . 
o  o 

Substituting  th'.s  in  Eq.  (l2)  gives 

Z  =  ap  exp[-h/K]  .  (l4) 

The  cutoff  frequency  for  an  ex*  onentially  tapered  transmission  line 
[Ref.  4]  is 


f 

co 


a 

4nH  ' 


(15) 


Since  both  the  velocity  of  sound  and  the  scale  height  vary  with  height 
in  a  similar  manner,  the  lowest  useful  frequency  is  not  strongly  height 
dependent.  An  average  \  ilue  of  0.004  cps,  taken  from  Eq.  (l5)  and 
atmospheric  characteristics  given  in  U.  S.  Standard  Atmosphere  [Ref.  3], 
was  selected  for  use  in  this  paper.  Figure  7  summarizes  the  atmospheric- 
attenuation  characteristics  for  a  sound  wave  traveling  vertically  from 
ground  level  to  220-km  altitude. 
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FIG.  7.  ATTENUATION  AS  A  FUNCTION  OF 
FREQUENCY . 

The  effect  of  atmospheric  filtering  on  the  particular  acoustic  wave 
f  interest — i.e.,  Brode's  6-percent  waveform — was  determined  by  analog 
means.  An  electrical  filter  having  the  atteiuation  characteristics  of 
Fig.  7  was  constructed.  When  excited  by  a  voltage  wave  shape  identical 
to  that  of  the  Brode  shock  wave  for  6-percent  overpressure,  the  output 
signal  represents  the  neutral  acoustic  disturbance  in  the  F  layer 
directly  above  the  blast.  The  total  acoustic-power  loss  factor  caused 
by  atmospheric  filtering--i ,e . ,  APf  of  Eq.  (l) — is  the  ratio  of 
maximum-output  to  maximum-input  voltage  (Fig.  8). 


APf  =  (V,/Va)2 


34273 

FIG.  8.  ANALOG  TECHNIQUE  USED  TO  DETERMINE 
POWER  LOSS  DUE  TO  SPECTRAL  FILTERING. 
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The  value  of  power-density  loss  factor  obtained  in  this  manner  was 


£P  =  4.75  X  10"3  .  (lfi) 

Thus,  when  the  values  from  Eqs.  (4)  and  (16)  are  substituted  into 
Eq.  (3),  the  power-density  ratio  is 

P2 

-  =  APg  x  ^ 

=  1.57  X  10  X  4.75  X  10" 

=  7.45  X  10"9  . 


The  value  of  F-layer  overpressure  may  now  be  determined  from  Eq.  (2) 
with  jt  =  6  percent  and  the  appropriate  values  of  y,  p,  and  a 
taken  from  U.  S.  Standard  Atmosphere.  The  predicted  maximum  F-layer 
overpressure  is  then  found  to  be 

«2  =  2.2  percent. 

The  above  calculation  has  ignored  the  effect  of  ground  reflections. 
For  the  case  of  perfect  reflection,  twice  as  much  energy  as  considered 
above  would  be  directed  upward.  Since  overpressure  is  proportional  to 
the  square  root  of  energy  density — Eq .  (l)  —  this  effect  could  increase 
the  predicted  F-layer  overpressure  to  a  value 

«2  =  3-1  percent. 

The  magnitude  and  wave  shape  of  the  predicted  ionospheric  pressure 
disturbance  are  shown  in  Fig.  9. 
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FIG.  9.  OVERPRESSURE  WAVEFORM  AT  220  KM  DETERMINED  FROM  ANALOG  MODEL. 

B.  INTERACTION  BETWEEN  ACOUSTIC  AND  RADIO  WAVES 

Acoustic  waves  traveling  in  the  lower  ionosphere  may  be  detected 
by  phase  measurements  on  hf  radio  signals  reflected  within  the 
disturbed  region.  Neutral-density  variations  associated  with  an  acoustic 

ft 

disturbance  result  in  corresponding  changes  in  electron  density.  The 
radio-wave  index  of  refraction  p  is,  of  course,  a  function  of  the 
local  electron  density  [Ref.  5];  the  phase  path  of  the  radio  signal  is 
thus  sensitive  to  any  neutral  pressure  changes.  The  rate  of  phase 
change  is  the  frequency  shift  or  "modulation"  observed  at  the  radio 
receiver . 

An  hf  radio  wave  vertically  incident  on  the  ionosphere  splits 
into  its  ordinary  (o)  and  extraordinary  (x)  components  (Fig.  10) 
as  a  consequence  of  the  geomagnetic  field.  The  raypaths  taken  by  the 
components  depend  on  the  frequency  used,  the  magnetic-field  geometry, 
and  the  structure  of  the  ionosphere.  Generally,  for  dip  angles  greater 
than  45  deg,  the  X  wave  deviates  less  from  the  normal  than  does  the 
0  wave . 

Both  components  are  sensitive  to  ionospheric  fluctuations  at  any 
point  along  their  paths.  The  greatest  sensitivity,  however,  occurs 

* 

In  the  absence  of  the  geomagnetic  field,  or  for  acoustic  waves 
traveling  in  a  direction  along  the  magnetic-field  lines. 
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FIG .  10 .  ORDINARY-EXTRAORDINARY 
RADIO-WAVE  SPLITTING  AT 
VERTICAL  INCIDENCE. 


near  the  point  of  reflection.  It  is  clea.  therefore,  that  the  exact 
phase  disturbance  resulting  from  an  acoustic  wave  propagating  through 
the  ionosphere  is  a  function  of  a  great  number  of  variables. 

For  the  purpose  of  calculation,  several  assumptions  were  used  to 
simplify  the  problem.  Only  the  X  component  was  considered,  and  ray 
bending  was  ignored  (Fig.  11 ).  The  ray  was  assumed  to  be  perpendicularly 
incident  on  a  parabolic  layer  being  perturbed  by  a  plane  acoustic  wave 
moving  at  an  angle  ¥  with  respect  to  the  vertical.  This  geometry 
simplifies  the  problem  greatly. 

The  electron  density  N  is  an  undisturbed  parabolic  layer  and  is 
described  by 


N 


(17) 


where  H  is  the  layer  half-thickness  and  h  is  measured  from  hQ 
(Fig.  12).  The  presence  of  a  pressure  wave  in  the  ionosphere  changes 
the  electron  density  by  an  amount  [Ref. 
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p-po  cos2  ♦ 
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where  <t»  is  the  angle  between  the  magnetic-field  lines  and  the  direction 
of  acoustic-wave  travel  The  resulting  new  density  profile  is 

N'  =  N  +  AN 


=-  N 

m 


/P-P 


2  \1 


1  + 


o  cos  n  ' 


7 


/J  ' 


(18) 


The  phase  path  of  a  vertically  reflected  radio  wave  is  found  from  the 
integrated  index  of  refraction 


h' 


(19) 


For  longitudinal  propagation  of  the  extraordinary  component,  the  index 
of  refraction  p  is  given  by  [Ref.  5] 


where 


and 


2 

|i  = 


1 


X 


„2 

‘N 

max  N ' 

2  N 
f  max 


f  cos  9 

n 

- 


(20) 


in  which  f^  =  maximum  plasma  frequency 
max 

fu  =  gyro  frequency 

n 

9  =  angle  between  magnetic  field  and  hf  path 
f  =  radio  frequency. 

Substituting  the  density  profile  of  Eq.  (18)  into  Eq.  (19)  yields  the 
expression  for  phase  path 
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where 

f(f  -  f  cos  9) 

C  =  - - -  .  (22) 

|l  +  [(p-P  )/p  ][cos2  0/7] 

N  I  O  O 

m 


The  upper  limit  h^  of  Eq .  (21 )  is  the  true  reflection  height  and  's 

found  by  setting  the  index  of  refraction  to  zc-ro — i.e.,  X  =  1  —  | Y  |  . 

L 

Or,  substituting  for  X  and  Y  ,  one  obtains 

L 


h1  =  h[  1  -  (l-C)/2l  . 


(23) 


Since  the  value  of  C  is  known  as  a  function  of  time  for  a  given 
radio-signal  blast -site  geometry  (Fig.  11 ),  the  radio-wave  phase  path, 
as  a  function  of  time,  may  be  found  by  substituting  Eq,  (23)  into  Eq . 
(2l),  The  total  phase  path  (3,  in  cycles,  is  then 


where  c  is  the  speed  of  light  and  f  is  the  radio-signal  frequency. 

Equations  (21 )  through  (24)  allow  one  to  compute  (3  as  a  function 
of  time  for  the  given  neutral  atmospheric  pressure  wave  (Fig.  9). 
Because  of  the  rather  complex  relationships  among  these  equations,  the 
computations  are  best  done  with  the  aid  of  a  digital  computer.  The 
approximations  involved  are  the  assumptions  of  l)  parabolic  ionospheric 
density,  2)  an  undeflected,  extraordinary  ray,  and  3)  a  pressure  dis¬ 
turbance  who  physical  size  is  much  larger  than  the  region  in  which 
the  index  of  refraction  differs  appreciably  from  unity — i.e,,  much 
larger  than  the  pressure-sensitive  portion  of  the  radio  path.  The 
equations  were  used  to  calculate  the  theoretical  phase-path  and  radio¬ 
frequency  disturbance  measured  at  Arneson  and  Bow  Island  during  the 
Snowball  test.  The  results  of  these  calculations  are  given  in  Chapter 
VI  of  this  paper. 
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III.  EXPERIMENTAL  EQU I PMENT 


The  ionosphere  above  the  Snowball  explosion  was  instrumented  with 
four  vertical-incidence  sounders.  Two  sounders  were  located  near 
Arneson,  Alberta,  85  km  magnetic  north  of  the  Suf field  Experimental 
Station;  the  remaining  two  were  installed  at  Bow  Island,  Alberta,  85  km 
magnetic  south  of  the  blast  (Fig.  l). 

To  obtain  the  greatest  possible  sensitivity  to  disturbances  in  the 
F  region  of  the  ionosphere,  the  pulse  transmitters  were  operated  in 
the  region  above  the  F-layer  ordinary-wave  critical  frequency.  These 
frequencies  were  measured  about  15  min  before  the  blast,  and  frequencies 
selected  at  that  time  were  used  throughout  the  test. 

Each  sounder  operated  with  a  peak-pulse  power  of  about  1  kw.  The 
pulsewidth  was  200  jisec ,  and  a  60-pps  repetition  rate  was  used.  Signals 
were  transmitted  and  received  on  separate,  horizontal,  three-wire, 
broadband,  dipole  antennas,  whose  centers  were  about  25  ft  above  ground 
level.  Figure  13  shows  a  plan  view  of  the  transmitting-  and  receiving- 
antenna  array  used  at  both  Bow  Island  and  Arneson. 


TRANSMITTING 

ANTENNA 

(HIGH  FREQUENCY) 


33302 


FIG.  13.  PLAN  VIEW  OF  TRANSMITTING  AND  RECEIVING  ANTENNA  ARRAY  USED 
AT  BOW  ISLAND  AND  ARNESON. 
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The  output  of  each  receiver  was  translated  down  to  25  kc  and  recorded 
on  one  channel  of  a  Sanborn  Model  2007,  seven-channel  magnetic  tape 
recorder  employing  standard  direct-record  electronics.  WWV  signals  were 
simultaneously  recorded  on  an  adjacent  channel  for  timing  purposes.  A 
block  diagram  of  the  system  is  shown  in  Fig.  14. 


34276 

FIG.  14.  BLOCK  DIAGRAM  OF  A  TRANSMITTER-RECEIVER 
PAIR  (two  at  each  site). 

Figure  15  shows  a  typical  signal  recorded  in  this  manner  for  a 
single  transmitted  pulse.  The  200-nsec  burst  of  25-kc  signal  derived 
from  the  transmitter  rf  pulse  is  followed  by  about  2  msec  of  delay 
corresponding  to  the  earth-ionosphere-earth  travel  time. 


SEL-65-062 


18 


25  kc  BEAT  FROM 
RECEIVED  ECHO 
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FIG.  15. 


A  TYPICAL  RECEIVER  OUTPUT  SIGNAL. 
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IV.  DATA  PROCESSING 


The  data  taken  at  Bow  Island  and  Arneson  were  recorded  on  magnetic 
tape  continuously  from  90  min  before  detonation  until  45  min  afterward. 
Figure  15  shows  the  recorded  signal  as  displayed  on  an  oscilloscope  for 
one  typical  transmitted  pulse  and  received  echo.  The  magnetic  tapes 
recorded  on-site  were  later  processed  at  Stanford  to  produce  range¬ 
time  film  records  showing  received  signal  phase. 

The  film  records  were  generated  in  the  following  manner.  An 
oscilloscope  sweep  was  triggered  using  the  leading  edge  of  the  25-kc 
signal  derived  from  the  transmitted  radio-frequency  pulse.  The  sweep 
was  then  intensity  modulated  with  the  remaining  recorded  waveform.  The 
result  was  a  line  across  the  face  of  the  oscilloscope  containing  a  series 
of  dots  corresponding  to  the  peaks  of  the  signal  shown  in  Fig.  15.  By 
moving  a  35-mm  film  strip  perpendicular  to  the  sweep,  the  dots  from 
each  successive  oscilloscope  trace  were  recorded  on  a  slightly  different 
portion  of  the  film.  The  film  was  moved  slowly  enough  to  allow  a  series 
of  pulses  recorded  in  this  manner  to  appear  as  lines  or  "stripes"  on 
the  film  strip. 

Figure  16  shows  a  typical  "stripe"  recording  taken  from  a  normal, 
undisturbed  ionosphere.  The  transmitted  pulses  may  be  seen  at  the  bottom 
of  the  recording  as  several  parallel  lines.  The  received  pulses  appear 
higher  up  on  the  record  as  diagonal  lines.  A  measure  of  the  phase 
change  introduced  by  ionospheric  fluctuations  may  be  obtained  by  counting 
the  number  of  stripes  that  cross  a  given  horizontal  line  through  the 
received  pulses.  Since  each  stripe  represents  a  2jt -radian  phase  change, 
the  relative  phase  in  the  path,  as  a  function  of  time,  may  be  plotted 
directly  from  these  records.  The  change  in  frequency  of  the  received 
radio  wave  is  then  simply  the  derivative  of  the  phase-time  record. 
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V .  RESULTS 


At  the  time  of  the  Snowball  test,  which  occurred  at  10:58  MST  on 
17  July,  1964,  all  four  pulse  sounders  were  operating  satisfactorily. 

The  two  Arneson  transmitters  were  placed  on  5.6  Me  and  5.7  Me  while 
those  at  Bow  Island  operated  on  5.2  Me  and  4.5  Me.  These  frequencies 
were  chosen  on  the  basis  of  the  measured  5.4-Mc  ordinary -wave  critical 
frequency  in  the  hope  that  the  return  echoes  would  contain  extraordinary- 
wave  components  only  (see  Chapter  II ). 

Clear,  explosion-induced  disturbances  were  detected  about  9'/s  min 
after  detonation  on  all  frequencies  except  5.2  Me.  Because  this  fre¬ 
quency  was  well  below  the  critical  frequency,  both  ordinary  and  extra¬ 
ordinary  waves  were  reflected  from  the  ionosphere  with  little  attenua¬ 
tion.  As  a  result,  the  return  echoes  from  successive  transmitted  pulses 
overlapped  and  produced  a  film-strip  record  in  which  the  F-layer 
perturbation  was  indi scernable . 

The  clearest  effect  was  measured  at  the  5.4-Mc  Bow  Island  site. 

Figure  17  shows  the  "stripe"  records  taken  immediately  before  and  during 
the  ionospheric  disturbance  on  this  frequency.  The  effect  is  character¬ 
ized  by  a  sudden  change  in  the  slope  of  the  return  echoes  and  is  most 
noticeable  at  about  10  min  after  detonation.  A  mixing  of  extraordinary- 
wave  (lower  return  echo)  and  ordinary-wave  (upper  return  echo)  signals 
occurs  at  about  10  min  30  sec.  Similar  disturbances  were  measured  at 
the  Arneson  site.  No  ordinary  components  appeared  on  5.6  and  5.7  Me, 
however,  because  the  transmitters  were  operating  about  200  kc  above 
0-wave  critical. 

The  phase-vs-time  record  taken  from  the  5.4-Mc  Bow  Island  "stripe" 
recording  is  shown  in  Fig.  18.  A  measure  of  the  effect  due  to  the 
explosion  was  obtained  by  drawing  on  the  phase  plot  a  curve  which  was 
considered  to  be  a  best  approximation  to  the  undisturbed  ionospheric 
return.  The  difference  between  the  two  curves  then  gave  the  phase 
disturbance  vs  time  record.  All  records  were  reduced  in  the  above 
manner.  The  results  are  summarized  in  Fig.  19. 
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FIG.  17.  STRIPE  RECORDING  TAKEN  AT  BOW  ISLAND  DURING  DISTURBED  PERIOD. 
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FIG.  18.  PHASE  VS  TIME  TAKEN  FROM  BOW  ISLAND  5.4-MC  STRIPE  RECORDING. 
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VI .  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 

The  purpose  ol  the  following  section  is  to  compare  the  experimental 
results  summarized  in  Chapter  V  with  curves  based  on  the  linear  acoustic 
theory  presented  in  Chapter  II .  Two  sets  of  on-site  data — the  atmospheric- 
wind  profile  and.  the  ionospheric-electron-density  profile--were  used  in 
calculating  the  theoretical  waveforms. 

The  estimated  horizontal  component  of  wind  blowing  in  the  Bow  Island- 
Arneson  direction  during  the  explosion  is  shown  in  Fig.  ?0.  Data  for 
this  curve  were  compiled  from  several  sources.  The  0-to-20-km  winds, 
for  example,  were  measured  with  a  weather  balloon  released  at  the  test 
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WIND  SPEED  (m/tec) 

34-280 

FIG.  20.  ESTIMATED  WIND  PROFILE 
OVER  SUFFIELD  AT  TIME  OF  BLAST. 
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site  shortly  before  T  .  Sounding  rockets  fired  at  Cold  Lake,  Aloerta 

o 

(about  300  miles  north  of  Suffield),  about  one  hour  after  the  blast  pro¬ 
vided  data  to  an  altitude  of  80  km.  No  measurements  above  this  height, 
however,  were  taken  during  the  test;  the  results  above  80  km  shown  in 
Jig.  20  were  estimated  on  the  basis  of  a  survey  of  several  articles 
giving  data  on  high-altitude  winds  [Refs.  7-10] . 

The  nearest  measured  electron-density  profile  at  the  time  oi'  the 
test  was  taken  at  Kenora ,  near  Winnipeg,  Manitoba,  with  an  NBS  C-4 
rounder.  Figure  21  shows  the  ionogram  from  this  station.  The  ordinary- 
wave  critical  frequency  is  seen  to  be  about  5.0  Me.  To  obtain  the  best 
estimate  of  electron  density  above  the  blast  site  at  Suffield,  Alberta, 
this  record  was  scaled  upward  in  frequency  so  that  the  critical  frequency 
corresponded  to  the  5.4-Mc  critical  frequency  measured  with  the  Stanford 
sounders  at  Bow  Island  and  Arneson  (Chapter  v). 

Radio-wave  and  acoustic  raytracings  (Figs.  22  and  23)  based  on  the 
experimental  wind  and  electron-density  profiles  outlined  above  were 
generated  using  digital  raytracing  capabilities  available  at  Stanford 
Electronics  Laboratories.  These  curves  were  superimposed  to  obtain  the 
onset  time  of  the  disturbance — the  time  at  which  the  acoustic  wave 
reaches  the  radio-wave  reflection  point — and  the  angles  <t  (the  angle 
between  magnetic  field  lines  and  direction  of  acoustic-wave  propagation) 
and  'i  (the  angle  between  the  acoustic-ray  direction  and  the  vertical). 

The  theoretical  radio-wave  frequency  disturbances  were  then  calculated 
with  the  appropriate  values  of  <J>  and  if  above  and  the  methods  given 
in  Chapter  II,  Section  B.  A  comparison  between  the  waveforms  obtained 
in  this  manner  and  the  experimental  results  is  shown  in  Fig.  24.  Table  1 
shows  the  onset-time  comparison. 
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FIG.  21.  VERTICAL-INCIDENCE  IONOGRAM  TAKEN  AT  KENORA ,  MANITOBA 


ALTITUDE  (km) 


34281 

FIG.  22.  RADIO-WAVE  RAYTRACING  BASED  ON  EXPERIMENTAL 
ELECTRON -DENSITY  PROFILE. 
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FIG.  23.  ACOUSTIC  RAYTRACING  OVER  SUFFIELD  BASED  ON  EXPERIMENTAL  WIND 
PROFILE . 


TABLE  1.  COMPARISON  OF  DISTURBANCE  ONSET  TIMES 


Field  Site 

Frequency 

(Me) 

Wave 

Experimental 
Onset  Time 
(Min ) 

Theoretical 

Onset  Time 
(Min) 

Bow  Island 

5.4 

0 

m 

+  9.5 

Bow  Island 

5.4 

X 

■ 

+  9.7 

Arneson 

5.6 

X 

+  9.5 

Arne son 

5.7 

X 

m 

+  9.8 
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VII .  CONCLUSIONS 

The  experiment  described  in  this  paper  has  demonstrated  the  effec¬ 
tiveness  of  using  a  relatively  simple  linear  acoustic  theory  to  describe 
the  propagation  of  explosion-induced  disturbances  to  large  distances. 
Figure  24  and  Table  1  summarize  the  comparison  between  results  obtained 
experimentally  and  those  calculated  with  the  linear  model.  The  agree¬ 
ment  is  quite  good,  particularly  in  Table  l--the  disturbance-onset-time 
results.  Inspection  of  Fig.  24  shows  that  the  period  of  the  predicted 
disturbance  is  slightly  shorter  than  that  measured  experimentally,  but 
the  wave  shape  and  the  amplitude  of  the  predicted  phase  disturbance 
compare  well  with  the  experimental  results. 

This  comparison  strengthens  the  confidence  in  the  sound-velocity 
profile  used  in  the  model  and  suggests  that  the  assumed  atmospheric- 
filter  characteristics  (Fig.  7)  may  be  in  error — that  is,  the  real 
atmosphere  appears  to  have  a  smaller  bandpass  region  than  shown  in 
Fig.  7.  Since  the  shape  of  the  filter  is  a  critical  function  of  the 
atmospheric  conditions  near  the  radio-wave  reflection  point,  it  may  be 
possible  to  obtain  some  information  about  high-altitude  constituents 
by  "fitting"  a  filter  to  match  the  observed  waveform. 

No  definitive  conclusion  about  the  effect  of  the  earth's  magnetic 
field  on  acoustic-wave  propagation  can  be  drawn  from  the  results  of 
this  experiment.  It  was  stated  in  Chapter  II  that  the  degree  to  which 
a  neutral  acoustic  wave  produces  a  corresponding  charged-particle 
disturbance  is  proportional  to  cos  $  (<J>  is  the  angle  between  the 

acoustic-wave  direction  and  the  magnetic  field).  In  an  attempt  to 
verify  this  dependence,  sounders  were  placed  85  km  magnetic  north  and 
south  of  the  blast.  Prior  to  the  experiment,  it  was  expected  that 
this  geometrical  arrangement  would  provide  a  $  of  about  30  deg  over 
the  south  site  and  60  deg  at  the  north  site.  Thus,  a  three-to-one  ratio 
in  measured  amplitude  was  expected. 

Unfavorable  northerly  winds  at  the  time  of  the  blast  (Fig.  20), 
however,  changed  these  angles  to  about  37  deg  and  41  deg  respectively. 
This  change  is  equivalent  to  a  l.l/l  amplitude  ratio  and  is  well  within 
experimental  error. 
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